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i/^ ■ Abstract At densities below the neutron drip threshold, a purely elastic solid model 

O , (including, if necessary, a frozen-in magnetic field) can provide an adequate description 

Q \ of a neutron star crust, but at higher densities it will be necessary to allow for the 

\ penetration of the solid lattice by an independently moving current of superfiuid neu- 

CT' trons. In order to do this, the previously available category of relativistic elasticity 
models is combined here with a separately developed category of relativistic superfiuid- 

• • . ity models in a unified treatment based on the use of an appropriate Lagrangian master 

• ^ \ function. As well as models of the purely variational kind, in which the vortices fiow 

^ \ freely with the fiuid, such a master function also provides a corresponding category of 

■ non-dissipative models in which the vortices are pinned to the solid structure. 
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1 Introduction 



The purpose of this work is to adapt the previously developed formalism for the treat- 
ment of various kinds of conductivity in a relativistic solid P to the particular case 
of a neutron superfluid in the solid crust of a neutron star. This treatment will in- 
clude allowance for the presence of a frozen-in magnetic field in the limit of perfect 
electrical conductivity, which will be a very good approximation in this context. (Such 
an M.H.D. type limit is more relevant in neutron stars than the kind of diamagnetic 
polarisation that was considered in earlier work The effects of magnetic fields 

may be significant for ordinary pulsars and will be particularly important in the special 
case of magnetars. In fact, the recently observerd quasi periodic oscillations in the 
aftermath of giant flares in the soft gamma ray repeaters SGR 1806-20 j3J |3] and SGR 
1900-1-14 jB] are suggested to be associated with oscillations originating in the crust 
after a major quake, see [3 El Ej and references therein. The magnetic field will be of 
crucial importance for a complete understanding of such oscillations. 

One of the main motivations for this work is to deal with the frequency glitches in 
pulsars, whose treatment has long been recognised to require allowance both for elastic 
solidity and for the independent superfluid motion of neutrons in the crust at densities 
above the "drip" threshold (at about lO^^g/cm^). In the literature on relativistic con- 
tinuum mechanics, the superfluid aspects (whose development is relatively recent) have 
mainly been treated separately [HI 1101 HH 1121 UH] from the purely elastic aspects, for 
which the formalism previously developed for that purpose [TH IT^ ITHl [T7j has already 
been applied [IHl UHl 1201 1211 1221 123 121| to fairly realistic kinds of neutron star models (as 
well as to idealised configurations [23 (201 of more artificial kinds that are of essentially 
mathematical rather than astrophysical interest). A synthesis of the kind that is needed 
j27| for the treatment of neutron conduction through the solid crust has however been 
recently developed in a Newtonian framework j2Hl; which should be accurate enough 
for the treatment at a local level of the outer and middle layers of a neutron star crust. 
Nevertheless, for the global treatment of a neutron star, and even for a local treatment 
of the deeper layers, a quantitatively precise description will require a fully relativistic 
treatment. 

As a step toward the achievement of this purpose, the present work will develop a 
treatment that is general relativistic in the sense of being applicable to a space time 
background with an unrestricted Lorentzian metric with signature { — ,+,+,+} having 
components g^j^i, with respect to arbitrary space-time coordinates x^. This metric will 
not be restricted by any specific evolution laws, so the treatment will be compatible not 
only with the standard Einstein theory of gravity but also with conceivable alternatives, 
and in particular with the trivial flat background case of special relativity, which will 
always be useful as an approximation in a local neighbourhood. We wish to stress that 
the formalism is not restricted to fixed backgound spacetimes. In particular, as an 
example, it is applicable to the full dynamical Einstein equations, for which it provides 
the right hand side. 

The Newtonian treatment referred to above j2Hl included allowance for dissipative 
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effects sucli as vortex creep, but - to avoid introducing too many complications at once 
- tlie relativistic treatment provided liere will be of conservative type, applying to the 
low temperature limit in which dissipative effects are neglected. As a generalisation of 
the variational formulation of the purely elastic relativistic solid model ^3] that is 
adequate for treating the outer crust layers (at densities below the "drip" threshold) 
the present treatment will be based on the use of an action density A of a kind that can 
be used as a Lagrangian to obtain models of strictly variational type, but that can also 
be used (in the manner to described in Subsection 13.21 ) as a master function for the 
construction of non dissipative models of more general type such as are capable, as has 
already been shown in the relevant multiconstituent fluid limit case JH] of incorporating 
the potentially important effect of vortex pinning. 

The category of conducting solid models developed here is designed for application 
on a macroscopic scale that is large compared with the intervortex separation. How- 
ever as a special zero vorticity case, apphcable on a mesoscopic scale - large compared 
with the microscopic scale of internuclear spacing, but small compared with the inter- 
vortex separation - this category includes models of the relatively simple hyperelastic 
subcategory [301 that is also of interest in a cosmological context. 

The requisite Lagrangian master function A will be formulated in terms of a solid 
structure, representing a distribution of atomic nuclei in terms of a 3-dimensional ma- 
terial base manifold with position coordinates g*, together with a current nf" of free 
baryons, representing superfiuid neutrons, and a current of confined baryons, repre- 
senting the protons and the fraction of the neutrons that are confined within the nuclei, 
so that their flow will be restricted by the comovement condition that is expressible 
(using a comma for partial derivation) by 

<q% = 0. (1) 



2 Relativistic action formulation 

2.1 Dynamic variables and background fields 

As well as depending on the "live" dynamical fields, consisting in this case of the 
(vectorial) current components n/^ and and the (scalar) base space coordinates 
and their derivatives, the Lagrangian will also depend on the metric g^i, which 
in the present treatment will be considered as a "dead" - meaning predetermined - 
background. In such a case the most general infinitesimal Eulerian (i.e. fixed point) 
variation of the Lagrangian will be decomposable as the sum of two contributions in 
the form 

5 A = (5^A + 5^ A (2) 

in which 6'^A is the realisable part attributable to a physically possible alteration 
of the configuration of the "live" dynamical fields, while 6^ A is a virtual part arising 
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from mathematically conceivable but (in the context under consideration) physically 
forbidden variations of the "dead" background fields that have been fixed in advance, 
of which the only one in the present instance is the fiat or curved spacetime metric g^iy, 
so that the complete background variation contribution will be given just by 

.*A^ (3) 

The particular kind of background variation to be considered in this section is one that 
is simply generated by the action of a displacement field C,'^ say, so that the ensuing 
field variations will be given just by the negatives of the corresponding Lie derivatives. 
For the background metric the resulting variation will be given (using round brackets 
indicate index symmetrisation) by an expression of the familiar form 

5g^, = -2V(^e.) , (4) 

in which denotes the (components of the) Riemannian covariant differentiation 
operator with respect to g^^. 

In so far as the "live" contribution is concerned, a further subdivision arises in cases 
such as those of interest here, which are characterised by a constrained variation prin- 
ciple, meaning one whereby the "on shell" evolution condition is that of invariance of 
the relevant action integral with respect to a compactly supported perturbation of the 
dynamical fields that is not entirely arbitrary but constrained by an appropriate admis- 
sibility condition. In the applications under consideration, the relevant admissibility 
condition will be interpretable as the requirement for the change to represent a "natu- 
ral" variation of the same given physical system, whereas more general changes would 
represent a replacement of the system by a slightly different system within the same 
category. In such a case the "live" field contribution will be decomposable (though not 
necessarily in a unique manner) as a sum of the form 

(5^A = 6^ A + (5«A (5) 

in which 6'^ A denotes a part the would be inadmissible for the purpose of application of 
the variation principle, whereas is a "natural" variation that would be allowed for 
this purpose. This means that for any (unperturbed) configuration that is "on shell" - 
in the sense of satisfying the dynamical evolution equations provided by the constrained 
variation principle - a generic admissible variation must satisfy the condition 

6^A = 0, (6) 

using the symbol = to indicate equivalence modulo the addition of a divergence (which, 
by Green's theorem, will give no contribution to the action integral from a variation 
that is compactly supported). 

In the kind of medium with which we are concerned, the constrained variables will 
consist of the free and confined current vectors n/* and n^'^, while the unconstrained fields 
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will consist of the triplet of scalar fields that are interpretable as local coordinates on 
the material base manifold, and from which, together with their gradient components 
q^^, an associated particle number (which might count ionic nuclei or crystalline solid 
lattice points in the relevant microscopic substructure) density current, rij say, can 
be constructed in the manner prescribed below. The corresponding "live" part of the 
generic action variation will thus be expressible in the form 

(9A 

S^A = ^i'^Snt + + VJq^ + g^Sq- , (7) 

in which f^i^ ^ and fi^ ^ are interpretable as the 4-momentum covectors of the free and 
confined baryons respectively. 

It is to be observed that the partial derivative V^^ = dk/dq^j^ has the notewor- 
thy property of being tensorial not just with respect to transformations of the space 
time coordinates but also with respect to transformations q^ i— > q^ of the material 
coordinates, q^ whose effect will simply be given by 

q-^r ^ (8) 

This good behaviour contrasts with the comportment of the remaining partial deriva- 
tive, dA/dq^, for which the corresponding transformation has the non-tensorial form 

gqA QqB QqA c ,v dq^dq^ 

For each of the separate currents and the admissible variations are gener- 
ated by flow world line displacements, which means that they will be specifiable by 
corresponding displacement vector fields, ^/^ and say, as 

d'^n^ = -£^^ - n^V, , = -^^c - C (10) 

in which the Lie derivative is just the commutator 

^6 < = ^f'^'^'^ - <v.ef^ = ^c^^< - (11) 

For the material base space, the corresponding the flow world line displacements will be 
given just by unrestricted variations of the base coordinates q^, which will be expressible 
in the form 

^y = -gtM^^ (12) 

in terms of their own displacement vector field, say, which can be chosen indepen- 
dently a priori, but which will subsequently be subjected to the identification = 
when the confinement restriction is imposed. 

It can be seen (using integration by parts) that, in terms of the displacement fields 
introduced by (jlO|) and (fT^ . the generic admissible variation for the current carrying 
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medium will be expressible - modulo a variationally irrelevant divergence term - in the 
form 

<5^A--/Uf'^-/U;-/'r, (13) 

in which the covectorial coefficient /^^ will be interpretable as the force density acting 
on the free neutron current, and f^^^ will be interpretable as the force density acting on 
the confined baryon current, while the extra coefficient is the supplementary force 
density due, as discussed below, to stratification or solid elasticity, that acts on the 
underlying ionic (crystalline or glass like) lattice structure. 

If we postulate that the system should obey the constrained variational principle 
that is expressible as the imposition of the "on shell" condition (jU)), then it is evident 
from (|13|) that the ensuing dynamical field equations will be expressible simply as the 
separate vanishing of each one of these force densities, i.e. as the requirement that we 
should have /^^ = 0, f^^ = and = 0. We shall however be mainly concerned with 
other possibilities, particularly cases in which one of the constituents, labelled "c" say, 
is "confined" in the sense of being subject to a convection condition of the form 
which means that it has to move with the underlying ionic flow. When the application 
of the variation principle is subject to the corresponding convective constraint 

Cc = C , (14) 

the ensuing system of dynamical field equations will no longer entail the separate 
vanishing of and f^^, but only of the amalgamated ionic force density f\ that is 
defined as their sum, 

/i = fl + fl , (15) 

which is interpretable as the net force density acting on the integrated system consisting 
of the ions in combination with the convected constituent. 

2.2 Noether identity for general covariance 

If we now replace the natural variations (fTIH) by the corresponding Lie variations 

Sn,^ = -£^^n,^, 5n^,=-£^X^ = (16) 

it can be seen that the current contributions will leave extra terms, so that we shall 
obtain a live variation of the form © with 5^K given by and with the extra 
"unnatural" (variationally inadmissible) contribution given by 

<5»A = + ^fj^ni; , 5^n^ = n^V.it , = n^V^^'^ . (17) 

We can now obtain a Trautman type Noether identity (that will hold regardless of 
whether or not the variational field equations are satisfied) by taking the variations 
of all the relevant fields, "live" as well as "dead", to be given by the action of the 
same displacement field C,'^ whose effect on the background was given by (@)), so that 
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the resulting effect on the Lagrangian scalar will be given just by the corresponding 
displacement variation 

SA = -CW^A ^ AV^r • (18) 
We thereby obtain a relation of the form 

r {fi + fl + fl) = + 5«A - AV.r • (19) 

in which, by Q, (|H) and (fTTj) . the "dead" and "unnatural" parts will be respectively 
given by 

^^A = -2t^ V,e. , = + f^innV.e ■ (20) 

The relation (jl9|) can thereby be rewritten in the standard form 

r {fu + n + fl) = -T^v.c = rv^n , (21) 

in which the geometric stress momentum energy density tensor can be read out as 

J.,. = 2 + <lfg^^^ , VP = A - - /.'^.n; . (22) 

Since the arbitrary field C,'^ can be taken to be non zero only in the immediate neigh- 
bourhood of any chosen point, the Noether identity (|2ip implies that, at each point, 
this stress energy tensor will satisfy a divergence identity of the simple form 

= fi + fl + fl ■ (23) 

In the "on shell" case for which the variational evolution equations are satisfied, the 
force density contributions on the right of (|23|l will all drop out. 



2.3 Canonical formulation 

If, instead of working out the variations in terms of force densities modulo a divergence 
as in (fT^ and ()21|) . we insert the displacement contributions fl20j) and (fTB|) into the "live" 
variation ((7j) directly in (j21) we get an identity in which the coefficients of the locally 
adjustable fields ^'^ and V^^'^ must vanish separately. The former of these conditions 
reduces to a triviality, but the latter provides a relation showing that the geometrically 
defined stress energy tensor ()23p can be rewritten in the equivalent canonical form 

n = + f^'X + f^%< - , (24) 

in which the stress contribution at the end is given in terms of the coefficient introduced 
just before (jH)) by 

V^ = V>{q%. (25) 



7 



The "live" variation formula (I7|) also provides what is needed for the derivation of the 
corresponding expressions for the force densities introduced in ()13|) which for the fluid 
constituents will have the form that is already familiar from the preceding work P . In 
terms of the corresponding generalised vorticity 2-forms defined, using square brackets 
for index antisymmetrisation, by 



w'^^ = 2V[^a/,] , w% = 2Vi,fil] , (26) 
these current force densities will be given by 

f. = <n^u + V^nf , n = + f^\V,n^^ , (27) 

while the force density acting on the underlying atomic structure of the medium will 
be given by the formula 

(28) 

in which the Eulerian derivative is given by a prescription of the usual form 

-V.V\, V\ = --. (29) 



§qA QqA A ' ^ dq 



A 
V 



It is important to notice that the non-tensorial base coordinate transformation property 
^ of the first term in this formula will be cancelled by that of the second, so that the 
effect on the Eulerian derivative of a change of the material base coordinates will be 
given simply by 

fA fA dq^ , , 

5gA SqB QqA 

which shows that the (stratification and solid elasticity) force density (j^Hj) is invariant 
with respect to such transformations, an important property that is not so obvious from 
its detailed expression 

fl = Q^t. + nvM^t. - V. (n'^tJ ' (31) 

in which only the final term is separately invariant. 

It is important to observe that, regardless of what - variational or other - dynamical 
field equations may be imposed, the force density (jHTj) will never be able to do any work 
on the medium since, as an obvious consequence of ()28|) it must automatically satisfy 
the further identity 

fx = , (32) 
where u'^ is the unit 4-velocity of the medium as specified by the defining conditions 

q^u" = , u^u"" = -c^ , (33) 
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where c is an adjustable parameter that is introduced to facihtate the derivation [221 
of the Newtonian hmit case. 

If the variation principle © were imposed for the three independent displacement 
fields, so that - as remarked above - by (fT^ the three separate force densities /^^, /"^ 
and would all have to vanish, then it evidently follows from the Noether identity 
that the divergence condition 

V^n = 0, (34) 

would have to hold. This condition is interpretable as an energy-momentum balance 
equation that must be satisfied whenever the system is not subject to any external 
forces. So long as we are concerned only with a system that is isolated in this sense, so 
that the energy-momentum balance condition (|34j) is satisfied, it follows from (j23p that, 
even for a more general model admitting the presence of dissipative or other internal 
forces, these must be be specified in such a way as to satisfy the total force balance 
condition 

fu + n + /I = , (35) 

in order for the model to be self consistent. More particularly, in view of ()32|) it can 
be seen that the fluid force densities by themselves will have to satisfy the condition 

«^(/Urj=0, (36) 
expressing energy balance in the local rest frame of the underlying medium. 



3 Conservative dynamical conditions 

3.1 Baryon conservation and the confinement constraint 

The system under consideration is designed for applications in which the total baryon 
current 

< = < + n^" , (37) 

can be decomposed as the combination of a free neutron current /if and a convected 
part nf consisting of the remaining baryons (some of the neutrons and all of the protons) 
that is convected with the nuclei, so that it is characterised by the confinement condition 

= ficu" , ric = {-ncun^f/"^ /c. (38) 

As a consequence of the condition of conservation of the total baryon current, 

V,n^ = , (39) 

it follows that in terms of the "free" current vorticity 2-form w^^^, the energy conser- 
vation identity (j36|) will reduce to the form 

u^u/^^ni = - V.n^ . (40) 
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Due to the confinement constraint, there is a freedom of adjustment in the specifi- 
cation of the confined particle momentum covector fi'^^ by the general formula ((Tj), but 
this does not affect the quantities involved in pUj) because there is no ambiguity in the 
specification of the energy component 

S^ = -u>^fi^^. (41) 

An obvious way of resolving the ambiguity in the remaining (space) components of 
the confined particle momentum would be to set them to zero, so that the complete 
covector would simply take the form /i^^ = —£^u^. Such a choice would be adequate 
for the specialised purpose matching the formalism used for the non-conducting solid 
limit case [23 • However for the broader objective of consistency with the formalism 
developed for the multiconstituent fluid limit case jU] we shall adopt a more satisfactory 
alternative ansatz for resolving the ambiguity in the specification of /ti'^^, which is to fix 
it by imposition of a convention to the effect that 

V\u, = Q. (42) 



3.2 Vortex pinning and chemical equilibrium 

As the condition of superfluidity of the "free" constituent will consist in a restriction of 
the possible behaviour of the corresponding momentum covector it is important to 
remark that this covector will not be affected by the chemical basis changes that will 
be discussed in Subsection I4.HI On a mesoscopic scale (large compared with the micro- 
scopic lattice spacing but small compared with intervortex separation) the condition 
of superfluidity will be expressed by the condition that this covector should have the 
form of a gradient, /i^^ = KV^ip of a periodic quantum phase angle ip, with the corollary 
that the corresponding "free" vorticity tensor u/^^ will also vanish at the mesoscopic 
level. At the macroscopic level, on a scale large compared with the separation between 
vortices, this vorticity will have a large scale average value that will not vanish, but to 
be compatible with a mesoscopic fibration by 2 dimensional vortex lines it must still be 
algebraically restricted by the degeneracy condition 

w'^^^w'^,^ = . (43) 



This condition will automatically be satisfied by the field equations that will result 
from the relevant variation principle, which will require that the vanishing perturbation 
condition (jUI) should be satisfied, not for independent variations of all three of the 
displacement fields ^Z*, which in view of the constraint ^7\i would lead to 

overdetermination, but just for variations satisfying the corresponding restraint ()14j). 
This requirement entails the vanishing, not of all three of the relevant force densities 
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/u5 f^fiy given by the prescriptions (j77j) and (jHH), but only of the free force density 
and of the amalgamated force density /^^ given by (fT3j) that acts on the part that 

is convected with the ionic lattice. 

According to the general formula ()27|) . the vanishing of the free force density /^^ 

entails not only the separate conservation of the free part n/* of the current density but 

also a dynamic equation of the familiar form 

= , (44) 

which is interpretable as meaning that the vortex lines are "frozen" into the free fluid 
in the sense of being dragged along with it. 

Instead of postulating the strict application of the variation principle, we can obtain 
a non-dissipative model of another physically useful kind by postulating that the vortices 
are "pinned" in the sense of being frozen, not into the fluid but into the atomic lattice 
structure, which means replacing by a dynamical equation of the analogous form 

<X = 0, (45) 

which unlike is unaffected by changes of the chemical basis [Oj that is used, as 
discussed in Subsection 14. 3t to specify which of the neutrons are to be considered to be 
"free" . 

Whichever of the alternative possibilities pij) or (jlSj) is adopted, the degeneracy 
condition will evidently be satisfied, and more particularly the left hand side of 
the identity pn|l will vanish, which means that the model must be such as to satisfy 
the condition 

- = . (46) 

In the strictly variational case this requirement will obviously be implemented by the 
separate conservation condition 

V,r7/ = 0, (47) 

which will be physically realistic in the context of high frequency oscillations. However 
in applications to slow long term variations it will be more realistic to use a model 
based on the alternative possibility, namely that of chemical equilibrium, as given by 
the relation 

- = , (48) 

which, like ()45j) is does not depend on which fraction of the neutrons is considered to 
be "free". 

Thus, depending on how we choose between the alternatives ()44j) or ()45|). and be- 
tween the alternatives ()47p or PHj) . we can use the same Lagrangian master function A 
for the specification of 4 different kinds of non-dissipative models, which are categoris- 
able as unpinned with separate conservation or chemical equilibrium, and pinned with 
separate conservation or chemical equilibrium. 
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4 Internal action function 



4.1 Subtraction of the ballistic contribution 

For the purpose of studying the Newtonian hmit it is convenient [211 to spht up the total 
action function into parts depending just on the confined number density ric defined by 
flHSj) and the corresponding free number density rif defined by setting 

< = nf< , n, = {-n,,n^Y/^/c , (49) 

in a decomposition of the form 

A = A,, + A,„. , (50) 

(using the tilde just to distinguish the parts that will be relatively large in the non 
relativistic limit) in which the first term is a ballistic contribution (meaning that by 
itself it would simply imply geodesic motion of the separate contributions) given in 
terms of a fixed mass parameter nio by an expression of the simple form 

A-bai = -"^o c^("'c + m) , (51) 

so that all the non-trivial information about the system will be contained in the re- 
maining internal action contribution A.^^, which remains finite in the large c limit [221 
(for which the tilde parts would diverge). 

As far as the relativistic treatment is concerned, mo may be given an arbitrarily 
chosen positive value or simply set to zero. The reason for taking the trouble of intro- 
ducing it at all is that in appropriate circumstances nio can be judiciously chosen in 
such a way as to make the internal remainder A.^^ relatively small, so that it can be 
dealt with as a small perturbation in an approximation of Newtonian type, for which 
mo will be interpreted as the "rest mass per baryon". In most astrophysical contexts 
for which a Newtonian approximation is applicable, and in particular for local pertur- 
bations in the crust of a neutron star, it will be good enough to use the proton mass 
calibration mo = nip, but it is well known that under laboratory conditions for which 
higher accuracy is needed it will be better to take the "rest mass per baryon" to be 
the standard atomic unit defined as a sixteenth of the mass of an ordinary O^^ oxygen 
atom. At the opposite extreme, at extremely high densities it might seem more natural 
to interpret the "rest mass per baryon" as meaning three times the "bare mass" of a 
quark, but subtraction of the relatively small rest mass contribution defined in that 
way would be pointless. In practise, under the highly relativistic conditions prevailing 
in the core of a neutron star, the contribution from A.^^^ will inevitably be relatively 
large, so that the most judicious choice may simply be to set mo = 0. 

However mo is chosen, its introduction will lead to a corresponding momentum 
decomposition of the form 
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in which the non-trivial "chemical" momentum contributions from the internal part of 
the action can be read out from its expansion in the form 



S'^K. = x'M + + + -Q^h' , (53) 

in which (since the specification of the ballistic part does not involve the q^) it can 
be checked that the coefficients in the last two terms will be the same as those already 
introduced in the analogous expansion ((7j) of the total action density. 



4.2 Material projections 

According to the general principle developed in the preceding work pP the internal 
action will be expressible, at a given material position, as characterised by the ma- 
terial coordinates g^, as a function just of the (time dependent) induced metric with 
contravariant components given by 

and of the relative current components 

= n.^qA = ^-^a ^ = ^ [^j. ^ ^^^u, = , (55) 

together with the corresponding material rest frame number densities, namely the 
confined density ric for which no projection is needed, and the free number density 
given in terms of the Lorentz factor 7 specified by the relative current velocity [v]'^ as 

mi = -c-\,n,' = 7 nf , 7 = (1 - M Vc')"'/' • (56) 



These conditions imply that the internal action density (which will be the same as 
the complete action density if we choose nio = 0) will have a generic variation of the 
form 

dA dA 

SA. , = -x^ Sn^ - Snn + p\ 5n^ + 7^^(57*^ + ^r^f^g* • (57) 

This provides an associated "convective" variation pQ (in which, with respect to ap- 
propriately dragged coordinates, both and are held constant) of the form 

dA dA 

8A= -x' Sn,-x^\Snii + —^5 nl + —^5 7^^^ , (58) 

c int A c C Ac"' Qjll^ Q'yfJ-'^ ell V / 

in terms of tensorial coefficients defined by 

dA. , dA. ,dA.^ 

ft7f=P^^ = PX- ifp-g^'AC- (59) 
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4.3 Ionic density and chemical gauge adjustments 

It is computationally convenient and physically natural to take the material base space 
to be endowed with a measure form that is specifiable in terms of antisymmetric com- 
ponents, ni^3(2 s^Y' ^^^^ fixed in the sense of depending only on the q^, and that 
will determine a corresponding scalar space-time field ni by the determinant formula 

ni = ^^^Iabc^Idef 7^^7^^7^^ . (60) 

For purposes of physical interpretation it will usually be convenient to take this measure 
to represent the number density of ionic nuclei, so that the confined baryon number 
density ric = —c~'^Uyn^ will be given by 

nc = Acni, (61) 

where Ac is the atomic number, meaning the number of confined baryons (protons plus 
confined neutrons) per nucleus. It is to be observed that, whereas the number current 
will fail to be conserved in the chemical equilibrium case characterised by ()48|1 . the 
formalism is such that the corresponding ionic number current 

= ni (62) 

must automatically satisfy the conservation law 

= (63) 

as a mathematical identity. Although there can be exceptions (for example when 
nuclear fission occurs) it will in fact most commonly be realistic to suppose that the 
number of ionic nuclei is indeed conserved even when the number of protons and neu- 
trons they contain is changing (whether by transfusive interchange of the confined part 
or just by conservative "dripping" of neutrons into the ambient superfiuid). 

For particular physical applications, the specification of the part of the baryon 
current that is to be treated as "confined" will be to some extent dependent on the 
timescales of the processes under consideration as compared with with the timescales 
for quantum tunnelling through the confining barriers containing the nuclei. For rela- 
tively rapid processes the fraction of the relevant quantum states that should be con- 
sidered to be confined may be subject to an increase, so the corresponding current, ric 
say, of baryons that are effectively confined will be somewhat larger than the value, nj^, 
that would be appropriate for processes occurring over longer timescales. It is therefore 
of interest to consider the effect of such a chemical basis adjustment, as given by a 
transformation of the form 

= (1 + e)nl , n'f" = < - en^ , (64) 

for some suitably specified dimensionless adjustment parameter e. So long as e is just 
a constant, it is immediately apparent, that such a transformation will leave the "free" 
4- momentum covector invariant, meaning that we shall have 

= fii , (65) 
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and, as in the fluid case [U] it is also easy to see that the entire stress energy tensor (j24p 
will be similarly invariant, 

n = n. (66) 



A less obvious observation is that this covariance property of the canonical 
stress energy tensor is not restricted to transformations for which e is constant. It 
will still hold whenever e is specified as a function just of the confined atomic number 
Ac = Tic/ni, and of the material position coordinates q^. Since it can be seen from 
()6m) that the variation of the ionic number will have a "live" (fixed background) part 
expressible in the form 

= m 7abO^''^<. + > (67) 

the corresponding variation of the adjustment parameter will have the form 

Unlike the free particle momentum which is subject to the chemical invariance condi- 
tion (jnSl), the confined particle momentum covector will undergo a non-trivial chemical 
adjustment given by 

/^^ = a; + - (e K - ^^f^^ ■ (69) 

However it can be verified that the extra terms will cancel out in the formula for the 
stress energy tensor so that the invariance condition will remain valid. 



4.4 Evaluation of momenta and stress energy tensor 

To proceed from (jHTjl . what we need is the corresponding "live" variation, as carried 
out at a fixed position in a fixed background, for which one can obtain an expression 
of the standard form 

6^A,^, = xl Sn: + xi + ^ I '^^f. + ^ Sq- , (70) 

in which, by using the constraint u^5q\, = —q^j^Su" the terms can be recombined in such 
a way as to obtain a well defined set of coefficients characterised by the requirement 
that should satisfy the condition (jl^ . It can be seen that the required (purely 
spacelike) value of P ^ will be given by 
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while the corresponding expressions for the internal contributions to the free and 
confined and particle 4-momentum covectors will be given by 

f X" II c ^ ^^'11 I I /^o^ 

Xu = ^u^+P Xu = ^u^ P u + ^—T^^u^ (72) 

which implies that we shall have 

X" = -Wx'u , = -u'xl ■ (73) 



It follows, according to the canonical formula ()24|). that the internal contribution to 
the stress energy tensor will be given by 

Tr..u = X'X + X^< + - V^u , (74) 

with 

n = n< + 2|^7^^ (75) 

and with the generalised pressure scalar 'J given, as in the multiconstituent fluid case, 
by 

* = A.. - X'X - X\< . (76) 
This can be rewritten in the manifestly symmetric standard form 

T^^ = p.^^^u^u^ + 2Q^i:y) + P^^: , (77) 

with 

Q>, = 0, Pr«. = 0, (78) 

in which the symmetric internal pressure tensor and internal momentum density vector 
will be given by 

n: = - 2 Y'r^^ , QL = x"nr , (79) 

while finally the internal mass density p.^^ will be given by an expression of the pseudo- 
Hamiltonian form 

cV,.=P".^r-A,.. (80) 
To put the total into the same standard form 

T^"' = pu^'u" + 2Q ^^u"^ + P'"' , (81) 

whenever the rest mass parameter ttIq is taken to be non-zero it will be necessary to 
use the correspondingly augmented quantities 

777 

PM- = ^n^^n1+ P[^^ , Q^ = m,rA< + <) + QL ' (§2) 
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so that, finally, we have 



p = rrio ( nc + — ) + p,^^ , (83) 



which can be rewritten as a Legendre type transformation formula giving the action 
density in terms of the rest mass energy density as 

A = {mo-f[v]^ + p\) nl - c^p . (84) 



5 Elastic energy contributions 
5.1 Relaxed contribution 

As in the multiconstituent fluid case [21] it is useful to decompose the internal action 
function in the form 

A. = A , + A , (85) 

mt cla ' ent \ / 

SO as to obtain a corresponding decomposition 

Pint = Pela + Pont , (86) 

in which A^^^^ and p^^^ are the parts attributable to entrainment ( which will be dealt 
with in Section lUI) and A , and p , are the static internal contributions that remain 

■—y cla ' cla 

when the relative current contributions are set to zero. For this static part (but not 
for the rest) the action density will just be the opposite of the elastic energy density € 
as defined by setting 

^ = P..c' , (87) 
which can be seen from (jHOl) to correspond to 

A, =-€. 

cla 



Whenever the comoving constituent can be considered to be a sufficiently good 
electrical conductor - as will be the case (on a macroscopic scale large compared with 
the scale of the relevant plasma frequency) in all the neutron star layers for which the 
present model is intended - it will be possible, in the manner recently described in 
the purely elastic case to include allowance for electromagnetic coupling just in 
terms of a "frozen in" field, with antisymmetric material base-space components, 
say, that (like the components rij^g^^ of the ionic measure) will be fixed functions of 
the material position coordinates q^. In such a case there will be a part, €!^^^ say, 
of the elastic energy density that depends on the not just directly but also via its 
dependence on the components F^b' which themselves depend on the q^, and it will 
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commonly be helpful to separate this magnetic part out for separate treatment (as 
described in Subsection 15. 3|) in a decomposition of the form 

^ = ^i„. + ^.a. , A.. = A.„. + A.a. , (89) 

in which the €!.^^^ is the strictly static remainder (the only part that would be needed 
if the material were an electric insulator) that will be left when the frozen components 
are set to zero. 

It will usually be convenient to further decompose the strictly static (meaning non 
magnetic) part of the elastic energy contribution in the form 

e'.. = g'o + ^.o, , = -Ao (90) 

in which & ^ is the part that remains in a relaxed configuration for which is min- 
imised for given values of the independent current components ric, rif, n^, and of the 
determinant of the induced metric 7^"^ . Fixing this determinant is equivalent to fixing 
the value of the conserved number density rii that is specified by (jUnj) . We shall use 
the notation 7^^ and 7^^ respectively for the corresponding relaxed values of 7^^ and 
its inverse 7^^ (as defined by 7ab7'''^ ~ ^a) which, for given values of nj, n^, rif, the 
energy minimisation occurs. Thus substitution of 7*'^ for 7^'' in the solidity term €^^^^ , 
or in the total G, will give a generically reduced value 

^.,<^.oM (91) 

but it will have no effect on the relaxed part, or on the ionic number density rij, for 
which we simply get 

€!r, = ^o^ ni = ni. (92) 



The relaxed contribution will evidently be of fluid type with a generic variation of 
the form 

= -6A^ = X II + x'^Sn, + X^Sn, - X^Jq- , (93) 

(in which the final term allows for the possibility of built in inhomogeneity in addition 
to the stratification due just to the variation of the atomic number ratio Ac) with 

Snj = ^nn^J^^^ + ^6q^. (94) 

The relaxed contribution to the pressure tensor as given by (|79|) will have an expression 
of the familiar isotropic form 

pAB^p AB^ p ,,,^^+^1 S S^J^^^C (95) 
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5.2 Solidity contribution 



Let us now consider the solidity contribution whose job is to allow for the effect (in 
^'j^^J of deviations of from its relaxed value 7^^^ (as determined by the scalars rif, 
ric, rii). Such deviations can conveniently be accounted for ^3] in terms of the constant 
volume shear tensor whose material base-space representation is specified as 

"^AB = i(7AB -7ab)' (96) 
which means that the corresponding space time tensor will be given by 



In most applications to behaviour of a perfectly elastic (rather than plastic or other 
more complicated) kind, it will be sufficient to use an ansatz of quasi-Hookean type 
|14j . meaning one in which the rigidity contribution has a homogeneously quadratic 
dependence on the deviation in the sense that it will be given by an expression of 
the form 

= i^^^'^^'^AB'^CD, (98) 

with 

2^ABCD _ 2^(AB){CD) — y^CDAB {^^) 

in which is the relevant solidity tensor, for which the check symbol is used 

to indicate that, for a given value of the material position coordinates g^, it depends 
only on the scalars rif and ric and (via ni) on the relaxed metric 7^^. The condition 
that 7^g and 7^^ must have the same determinant entails that on the 3 dimensional 
material base the symmetric shear tensor will have only 5 (instead of 6) independent 
components, and more specifically that to first order it will be trace free with respect 
to either the actual metric 7^^ or the relaxed 7^^^. It is therefore necessary to impose a 
corresponding restriction to completely fix the specification of the solidity tensor E^^'^^ , 
which can be conveniently done ^3] by requiring that it be trace free with respect to 
the relaxed metric 



It is to be remarked that an alternative procedure that would be equivalent for 
most practical purposes would be to work in an analogous manner with a shear tensor 
constructed as a difference between contravariant rather than covariant versions of the 
actual and relaxed forms of the induced metric. Another similarly equivalent possibility 
would be to use a judicious compromise between these alternatives, of the kind that 
has recently been shown jSHI to be technically advantageous for certain purposes. In 
the small shear limit that is physically relevant, the differences between such alterna- 
tives are only of quadratic order 0{(;^} for the shear itself and of higher order for the 
energy density €!^^^ , and will therefore usually be too small to matter in practise, so the 
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question of whether to use a prescription of the originally proposed kind ^1] or of a 
more sophisticated alternative kind |221 can be made just on the basis of mathematical 
convenience depending on the context of application. 

The specification of the solidity tensor is not by itself sufficient to complete the 
specification of the elastic system, as it is also necessary to specify the dependence on 
rii of the relaxed inverse metric 7^^^. The simplest possibility is that of what has been 
termed a perfect solid ^3], meaning one for which the elastic structure at each material 
position is isotropic with respect to the relaxed metric, which in that case can vary only 
by a conformal factor. This means that it will be given in terms of its value 7,^^ say at 
some fixed reference value no say of the ionic number density rij by 

^AB = (n,/no)'/'7o^" , 7^B = (^o/ni)'/=^7oAB (101) 



In the case of a solid structure that is isotropic (as will typically be the case on 
a macroscopic scale after averaging over randomly oriented mesoscopic crystals) the 
solidity tensor in the quasi Hookean ansatz will simply have to be given in terms of the 
relevant scalar shear modulus E by the formula 

^ABCD ^ 2IJ ('jMC^O)B _ l^AB^CD^ ^ (^^02) 

SO that ()98p will give the simple formula 

A., = -€^., = -i:?^ (103) 
in which the scalar shear magnitude <^ is defined by the formula 

in which the final term will in practise be negligible since of quartic order, 0{(;^}, in the 
small <^ limit that is relevant, as the trace is already of quadratic order, 7^^'^ab ~ ^{"J^}- 
Under these conditions the solidity contribution will provide a pressure tensor given by 
the formula 

pAB = _^ABCD,^^+p^^^AB^ (105) 

of which the final term is a pressure contribution given by 

Pso, = nt— + n,— + m— + -]<,\ (106) 



(9nf dric drii 3 

which, since it is of quadratic order in <;^, will be negligible compared with the first 
(linear order) term for most practical purposes. (It is to be noted that the sign of the 
last of these quadratic order terms was one of the errata in the original treatment 
^3].) The corresponding pressure adjustment contribution for the canonical formula 
(|7^ will be given by 

= 7.a1;^--.p. - + f^) , (107) 

in which, again, the quadratic order term at the end will be negligible in practise. 
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5.3 Magnetic contribution 

The presence of a frozen in magnetic field (whose effect is hkely to be significant in 
apphcations to ordinary pulsars, and of paramount importance in the particular case 
of magnetars) is representable [221 by the specification of a fixed and closed material 
base-space two-form field with components -Fab- The closure condition means that it 
will be locally expressible in terms of a fixed gauge 1-form field with components 
on the material base manifold as its exterior derivative 

^AB = 2A[B,A, (108) 

The corresponding space-time tensor field 

F,u = F,^q-ql (109) 

will therefore be given by 

F^, = , A, = A^qf^ (110) 

and will automatically satisfy the perfect conductivity condition 

F^X = 0, (111) 

This condition is interpretable as meaning that field has no electric part, and so will 
be entirely determined by its magnetic part which will have material base components 
given it terms of the alternating tensor e^^'^ associated with the induced metric 7^"^ 

by 

B- = y-^F^^. (112) 

If polarisation effects can be neglected, which will probably be a good approxima- 
tion in the crust (though not the core) of a neutron star, then the magnetic energy 
contribution will just be the negative of the usual vacuum action density associated 
with the field F^^, which in terms of its materially projected components gives 

^.a. = , B' = i,,B-B- = i7^-7^°i^Ac^Bn • (113) 
from which one obtains 

= — (fi^fiB - fiV^) • (114) 

mag ^yj- ^ ' ' ^ ' 

It follows that the magnetic pressure tensor will be given by 

^ntg = ^(^V^ - 2B^B^) , (115) 
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so that in terms of the magnetic field covector 

B, = -f,,B-q-^, (116) 
one obtains an expression of the famihar form, 

T^"" = — ( ^uf'u'' + B^Y'^ - 2B^'B''\ , (117) 

for the magnetic stress energy tensor. 

6 Effect of the current 
6.1 Entrainment contribution 

In general the entrainment action function A^^^^ will depend, for a given values of q^, 
on the relative current components = g^rix as well as on the scalar magnitudes rif , 
Tic, and the induced metric components 7^^ so its generic variation will have the form 

SK.. = - Xlfco + 1^57" + %^Sn^ + ?^5,^ . (118) 

The entrainment action A , is characterised by the condition that it vanishes when 

ent ^ 

the relative current components are set to zero, so when these components are 
sufficiently small, as will typically be the case, it will be a good approximation to take 
this contribution to have the homogeneous quadratic form 

A , = — n^n^ , (119) 

272 II AB ' ^ ' 

in which the entrainment mass tensor has components rn^^ that, like the static action 
contribution, are independent of the current components n^, so that the corresponding 
partial derivative in pi8j) will be given by 

OA , 1 ± , , 



It is conceivable that the relaxed action function might involve a built in anisotropy 
favouring relative currents in some particular direction, but in cases of the simplest 
kind, to which the remainder of this subsection and the next will be restricted, this 
function A^^^^ will be of purely fiuid type in the sense that for given values of it will 
depend only on the set of four scalar magnitudes consisting of rii, rif together with 
the relative current magnitude r2± that is defined in terms of the (unrelaxed) metric 
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value which - using the material index lowering operation specified by the induced 
metric 7^^ - will be given by 

Such a functional dependence provides an expansion 

-^Ae. = -xlM - Xljn. - xlJni + X^q^ + ^5nl , (122) 

of similar form to the perfect fiuid contribution (j93|l but with an extra term involving 
a partial derivative that provides an expression of the form p2()|l in terms an isotropic 
mass tensor given by 



I r r dA 



while the other partial derivatives in ()118|) will be given by 

dKnt 1 f ,0 \ aA,„, . ^ dni 



f , s ^ "^A.^, g g dni 

^A^B +X.^i7.B J , ^ = -X.^- (124) 



The scalar m^^ introduced in this way is identifiable as the increment m^^ = — nio of 
the effective mass of the free baryons (meaning the superfiuid neutrons) as compared 
with the ordinary baryonic mass ttIq. This mass increment is expected to be positive 
(and in some layers large) [HI] in the solid neutron star crust, but (moderately) negative 
in the fiuid layers below. 



6.2 Relaxed action contribution 

It will be useful for many purposes to replace the decomposition (j85|) of the internal 
action density by an alternative decomposition of the form 

A. =A, +A., A.=A,+A , (125) 

int lax rig ' rig sol ' mag ' V / 

in which the relaxed - meaning shear independent - part will evidently consist of the 
combination 

A,. = Aq + A,„, . (126) 



The use of such a combination is particularly convenient whenever the entrainment 
contribution is of the isotropic type characterised by the variation expansion ()122j) . 
in which case the complete relaxed action density will have a variation given by an 
expansion of the analogous form 

f 

^K. = -xlM - xL^-c - xlM + ^LM + ^^^'^ (127) 
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with 

= +Y^ Y*" = y'' ^y'' Y^ = Y^ +Y^ = + A^ fl28) 

^lax '^'O ' ^lax '^O '^<=nt ' '^lax '^O '^<=nt ' lax A QA ' out A " \^^^} 

It can be seen that the relaxed contribution to the "hve" action variation (j7(J|) will 
simplify to provide an expression of the form 

=x! u + u ^< - x^(5^ni + A s , (129) 

lax '\z.yiV t -^laxl' C /V ^ ' lax A ^ ' ^ ' 

in which u ^"^^ v be read out as the relaxed parts of the internal momenta 
given by (ff^ . Since we can write (?b'^n\ = —u^S'^n^, it can be seen that, as the 
analogue of these free and confined particle 4-momentum contributions, we shall also 
be able to read out a corresponding ionic stratification 4-momentum contribution y, 
so as to obtain a complete set of relaxed internal momentum covectors that will be 
expressible in terms of the relative flow velocity vector introduced in ()55|) by 

= ^ ^'^ + '^^ c[v]y . (130) 

X.aV = + ^ - M.^ (131) 

XL^ = , (132) 

It follows that the corresponding contribution to the generalised pressure scalar (|76|l 
wiU be given by 

^lax = A,, + + ^cX,l - "^cM' • (133) 

In the formula (jTlj) for the extra stress, it transpires that the contributions involving 
the mass increment m^^ (proportional to dA^^^/dnl) will cancel out, leaving just the 
part due to stratification which will be given simply by 

K., = -xL^i^^- (134) 

It is to be observed that the need to include this extra term in the canonical stress 
energy formula can be avoided by treating the stratification momentum (jl32p as if it 
were on the same footing as the others, which means including a corresponding extra 
term in a supplemented pressure scalar defined by 

= Aax - - - X^nr = + x,!n, . (135) 

Combining the relaxed internal contribution A, with the ballistic contribution A. , 

O lax bal 

one obtains the complete liquid (meaning hydrodynamic as opposed to solid or magne- 
tohydrodynamic) part of the action density in the form 

A,. = A^ , + A, = A - A . , (136) 

liq bal ' lax rig ' V / 
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for which the corresponding hquid contributions to the momentum covectors (j52p will 
be given (in terms of the relative flow velocity [vY = n'^^n'[) by 

In terms of these quantities the corresponding liquid contribution to the stress energy 
tensor will be given simply by 

K'^ = ^'LX + + m'.^c^ + ' (138) 

and the corresponding elastically relaxed force contributions acting on free and confined 
currents will be given by an ansatz of the standard form (j27|) which gives 

4, = 2nf V[^a4.,.] + A^.V^n/^ , (139) 

and 

= 2n,'^V[^//;,] + ^^^.W,n^ , (140) 

while for the analogously defined force contribution (j2Hl) due to stratification acting on 
the underlying ionic lattice, it can be seen that similar reasoning leads to an expression 
of the slightly different form 

f = 2nf Vr. .1 + A , A ^ = A ^ q\. (141) 



As remarked at the outset, because the confined particles and the ionic medium are 
constrained to move together the confined and supplementary force densities are not 
separately meaningful: all that is relevant for the equations of motion is their sum, 
the amalgamated ionic force density (|15p. Its liquid contribution can conveniently be 
expressed by introducing the corresponding amalgamated ionic 4-momentum that is 
definable as 

/^uq- = /f^iiq.' + Ac/t4., (142) 

(where = Uc/ni is the atomic number as introduced above). In terms of the 
amalgamated ionic 4-momentum it can be seen that one will simply obtain the formula 

= 2"i"V[.mI.] + C.V.A, + , (143) 

while the associated contribution to the stress energy tensor ()24|1 will be expressible 
by the neatly concise formula 

Tr^u = Kf^l. + + , (144) 



in which we shall have 



lax t laxP i laxP ^ 
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6.3 Complete description for perfect conducting solid. 

Let us conclude by putting all the pieces together for the case of what we refer to as 
a perfect conducting solid - meaning a model of the quahtatively simplest kind whose 
structure is fully isotropic with respect to the relaxed metric. In so far as the liquid part 
is concerned, this condition means not only that the incremental mass matrix has to be 
isotropic, as supposed in the preceding subsection, but also that the extra stratification 
vector in ()141|) must vanish, A ^ = 0. 

Replacing the relaxed contribution ()142p in ()144p by the corresponding total ionic 
4-momentum covector given (in terms of the confined atomic number Ac) by 

= 1^1 + A,i2l , (146) 

we can apply an analogous tidying up operation to the complete stress energy tensor 
fl24j) which will thereby acquire the form 

T^, = (A - - n,'f^\)6'^, + nl'^i + n,y\ + T^^^, - V^^^ , , (147) 

in which the magnetic contribution is given by ()117j) and the only manifest allowance 
for the effects of solid rigidity is in the final term. 

If we make the further realistic (and chemically invariant) supposition that (since it 
arises just from Coulomb interactions between the confined protons) the shear modulus 
depends only on ric and nj, then the solidity term will not contribute to the free particle 
momentum so we shall have 

= , /i = fl. , (148) 

while the force density acting on the part that is comoving with the ions will be given 
by 

f\ = flu + fiu , (149) 

where the rigidity term (allowing for the combined effect of ionic solidity and the 
magnetic field that may be frozen into the accompanying degenerate electron plasma) 
will be given by 

/I = V.T^^ , T^=T''^ + T'^^ . (150) 

J rigf rigf ' rig!^ magi' gol'' ^ ' 

To linear order in the shear amplitude - which will be sufficiently accurate for most 
practical purposes - the solidity contribution will be given by the expression 

T:. = -V:. +0{^}, (151) 

with ^, given by the perfect solid formula p07|) . which will itself be expressible to 
first order in the shear amplitude <^ by a prescription of the familiar form 

V>:^, = 2IJ,^, + 0{e}, (152) 

in which IJ is the ordinary shear modulus (which elsewhere is commonly denoted by 
the symbol fi that, in the present context, has already been used for the designation 
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of the material momentum components). This means that for practical purposes, for a 
conducting solid of the perfect (meaning intrinsically isotropic) type the term that is 
needed j2Zj to allow for deviation from behaviour of (multiconstituent) fluid type will 
be entirely contained in the magnetic contribution of the usual magnetohydrodynamic 
form given by (|117|) and the extra solidity term as given by (|152j) . on the understanding 
that S is prescribed function just of Uc and rij. 
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